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Abstract 

An analytical nonlinear description of field-line wandering in partially statistically magnetic 
systems was proposed recently. In this article we investigate the influence of the wave-spectrum 
in the energy-range onto field line random walk by applying this formulation. It is demonstrated 
that in all considered cases we clearly obtain a superdiffusive behaviour of the field-lines. If the 
energy-range spectral index exceeds unity a free-streaming behaviour of the field-lines can be found 
for all relevant length-scales of turbulence. Since the superdiffusive results obtained for the slab 
model are exact, it seems that super diffusion is the normal behavior of field line wandering. 

PACS numbers: 47.27.tb, 96.50.Ci, 96.50.Bh 
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I. INTRODUCTION 



Understanding turbulence is an issue of major importance in space physics and astro- 
physics; see, e.g., in Refs. It has been demonstrated in several articles 
that stochastic wandering of ma gne tic field-lines directly influences the transport of charged 



cosmic rays (see e.g. BIB ,0,111 
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13|,Q,[l5|. Several theories have been developed to 



describe field-line random- walk (FLRW) analytically. The classic work of Jokipii (see 16|). 
for instance, employed a quasilinear approach for FLRW. In this theory the unperturbed 
field-lines are used to describe field-line wandering by using a perturbation method. It has 
often been stated that this approach is correct in the limit of weak turbulence where it is 
assumed that the turbulent fields are much weaker than the uniform mean field (SBi <C B ). 
To achieve a more reliable and general description of field-line wandering Matthaeus et al. 
(see [l7]) developed a nonperturbative statistical approach by combining certain assump- 
tions about the properties of the field-lines (e.g. Gaussian statistics) with a diffusion model. 
More precisely, in the Matthaeus et al. theory of field-line wandering is has explicitly been 
assumed that field-line wandering behaves diffusively. 

An improved theory for FLRW, which is essentially a generalization of the theory of 
Matthaeus et al, was recently developed by Shalchi & Kourakis (see 18J). By explicitly 
assuming a diffusive behavior of the field-lines, the Matthaeus et al. theory can be obtained 



from the S 
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chi & Kourakis approach as a special limit. However, it has also been demon- 
strated in [18|] that for slab/2D turbulence geometry, the field-lines behave superdiffusively. 
Thus, the Matthaeus et al. theory cannot be applied for slab/2D composite geometry. As 
also demonstrated in [3], quasilinear theory is only correct for pure slab geometry or for 
small length scales. 

In most past studies a constant spectrum in the energy-range has been assumed (in this 
case the energy-range spectral index is equal to zero). It is the purpose of this article to 
explore different values of the energy-range spectral index. The layout of this article goes as 
follows. In Section 2, we discuss different forms of the wave-spectrum which are appropriate 
for solar wind turbulence. In Section 3, we calculate the FLRW for pure slab geometry for 
different values of the energy-range spectral index by applying the exact formulation for 
field-line wandering. In Section 4, we employ the nonlinear theory of Shalchi Sz Kourakis for 
FLRW, in order to deduce an analytic form for the field-line MSD for pure 2D turbulence. 
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These results can easily be combined with the pure 2D result to describe field-line wandering 
in the slab/2D composite model (Section 5). In Section 6 we summerize our new results. 



II. DIFFERENT FORMS OF THE WAVE-SPECTRUM 



In [19| a two-component turbulence model has been proposed as a realistic model for 
solar wind turbulence. In this model we describe the turbulent fields as a superposition of 
a slab model (k || Bq) and a 2D model (k _L Bo). In this case the xx— component of the 
correlation tensor can be written as 



puk) = p s x l :\k)+p^{k) 



2D I 



(1) 



with the slab contribution 



and the 2D contribution 



yslab 



(k) = g**(k\\) 



Ah) 



In previous studies the forms 



gS i ab{h) = 9R L]nh5B 2 



2vr 



L s labOr> slab \ 



1-S 

k 2 



k 2 ! 2 ) 



for the slab wave-spectrum, and 



g 2 D{k±) = 2 _£M 12D sb 2 2D (i + k\i\ D r 
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(2) 



(3) 



(4) 



(5) 



for the 2D wave-spectrum were used. The energy-range of the spectrum is defined for 
fey I sil b and k±_ Z^. Clearly, both spectrum forms are constant in the energy-range. 
However, as discussed in several previous articles (see e.g. [3]), we find in heliospheric 
observations a steeper spectrum (according to [3] the energy-range spectral index - cf. (jSJ) 
below - should be q = 1.07). In the following we deduce and discuss analytical forms of the 
wave spectrum for slab and 2D turbulence models. 
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case Normalization constants Cj 
0<q<l Cl : =( T l 5 + s i T )" 1 

« = 1 c 2 : =i(K^ri^)) 

— 1 / \ qx 

1 < g < 2 c 3 := L ^fc min Z sZab J 

TABLE I: The exact values of the various normalization constants c% (i = 1,2,3) are provided. 
These expressions are correct for k m i n l s i a b <C 1. 



A. General form of the slab wavespectrum 

According to solar wind observations, the following form of the spectrum should be 
appropriate: 



g^ih) = ^slabSB 



2tt 



slab 



X < 



if 

I slab 

(k\\ l slab )~ 2v if 



h < k 



9 if U . < U.. < 7-1 
11 ftmin ^ I' slab 



I slab < % 



(6) 



Here we defined the slab-bendover-scale l s i ao , the strength of the turbulent field SB 2 lab , and 
the inertial-range spectral index 2v. k m i n indicates the smallest wave-number which might 
be related to the bulk plasma length scale L via k m i n ~ L~ x . We have also introduced the 
energy-range spectral index q. By taking into account the normalization condition 



5B 2 slab = 5Bl + 8Bl 



d 3 k 



p: i : b (k)+p^ b (k) 



(7) 



we find for the normalization constant Cj the values shown in Table [B The values shown 
there are valid if the condition k min l s i ao <C 1 is fulfilled. 
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case 


Normalization constants di 


< q < 1 
q = l 
Kq<2 


d-2 ■— ^ n (fc min z 2D ) 
d 3 := (q - l)(k min 


-1 

r 

fee) 



TABLE II: The exact values of the various normalization constants di {% = 1,2,3) are provided. 
These expressions are correct for k m i n l2D <C 1. 



B. General form of the 2D wavespectrum 



For the 2D spectrum we can adopt the same form for the spectrum as used in the last 
subsection for the slab spectrum: 

di 



g 2D {k±) = ^ kD5B 2 



2tt 



2/> 



X < 



if k±< k min 

(k± hD)~ q if fcmin < k ± < ZjjJ 



(8) 



(k ± l 



2D) 



-2v 



if 



1 2D < k±- 



Here we used the 2D-bendover-scale Iidi the strength of the turbulent field SB^d, and the 
inertial-range spectral index 2v. k m i n indicates again the smallest wave-number and q is 
again the energy-range spectral index. Fulfilling the normalization condition 



5Bl D = SBl + SBl 



d 3 k 



PL D (k)+P™(k) 



(9) 



we find for the normalization constant di the values shown in table [TT1 In the following 
we consider different values of the energy-range spectral index q and calculate the field-line 
mean square deviation analytically for pure-slab, pure-2D, and two-component turbulence. 



III. FLRW FOR PURE-SLAB TURBULENCE 



As shown in several previous papers (e.g. [181 ]) the field-line mean square deviation can be 
calculated exactly for pure slab geometry. For standard forms of the wave spectrum, where 
q = 0, we find the classical diffusive result: ((Ax) 2 ) = 2kft,\z\ (see e.g. [if], Q), with 
the field-line diffusion coefficient kfl- Shalchi & Kourakis (see 18[]) derived the following 
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ordinary differential equation (ODE) for the mean square deviation and slab geometry 



dz 1 



;(Ax) 2 > 



d A kP s x i :\k) C os(k\\z) 



B 2 

o_ P<X> 

^ / dk\\ g sla \k\\) cos(k\\z) 



B, 



(10) 



JO 



For the wave spectrum of Eq. ([6]) this becomes 

J2 r D 2 rl7,\ 



d 2 n , <)B 2 

± <(Ax) 2 ) « 4c^^§M / d*,,^ 

*^ k-min 



' cos(k\\z) + . . . 



(11) 



It can easily be proven that the contribution of the inert ial-range (k\\ > lJi a b) is much smaller 
than the contribution of the energy-range (k\\ < Z^ ob ), and was thus neglected in the right- 
hand side (rhs) of Eq. fTTTl) . Furthermore, the upper limit of the fen-integral can be extended 
to infinity. Here q > is assumed, for convergence. 
Taking into account the relation 



noo i 

/ dx x^ 1 cos x = — [ e - m ^ 2 T (/i, +iu) + e+^ /2 r (/i, 

Ju ^ 



-IU 



(12) 



according to Gradshteyn & Ryzhik (see [211 ]. page 430, Eq. 3.761.7 therein), for // < 1 
(implying here q > 0), where we have employed the incomplete Gamma function r(/i, x) = 
dt t M_1 e - ' (see Eq. 8.35 in the latter reference), and approximating r(/i, x), for small 
values of the argument 



v(n,x < i) ~ r(/i) 

(see Eq. 8.354.2 in the same reference) we find 



x h 



f°° f 71 \ 1 

/ dx x^' 1 cos x ~ T (/i) cos ( — ji j 



(13) 



(14) 



By applying this formula onto Eq. (TTTj) one gets 
dz 2 



;(a^) 2 > 



A R 2 

, 1 l-g ur> slab ^q-1 



^slab 



Bl 



r (1 - g) sin (y) + -^—^{zk r , 
The result can easily be integrated to obtain 
<(Ax) 2 > 



(15) 



4c, 



jl-q dB 2 lab +1 

<z(<Z + 1) s/ab £ 2 



1) 



(16) 



This expression is valid for < q < 1 and for 1 < q < 2. For q = 1, Eq. (TTTj) can be directly 
evaluated and we find a logarithmic behavior of the MSD. In the following, we shall further 
simplify Eq. ( TTB1) . by distinguishing the ranges < q < 1 and 1 < q < 2. We stress that the 
are interested in the large z range, although we note that the condition zk m i n = e 1 is 
assumed to hold everywhere (since k min is related to the inverse size of the plasma "box"). 
We therefore retain the definition of the small parameter e, whose polynomial contribution 
may be singled out, for order of magnitude estimates. 



A. Smooth spectrum form: the case < q < 1 

In this case the first term in the rhs of Eq. ffl6|) is dominant and we obtain 

«^ 2 > * wh 1 -^ 1 ^ (1 - »> sm (t) ~ • (17) 

while a contribution ~ 0(t l ~ q ) within the brackets in (fl6j) is omitted. In general the mean 
square deviation of the field-lines has the form {|(Ax) 2 ) = az b . According to Eq. (ITT)) we 
find for the slab model and for the values of the energy-range spectral index considered the 
characteristic exponent 

b = q + l. (18) 
It is obvious that we obtain super diffusion (1 < b < 2) for < q < 1. 



B. Steep spectrum form: the case 1 < q < 2 



In this case the second term in the rhs of Eq. (fl6l) is dominant (of the order e 1 q 3> 1) 
and one gets 



;(Ax) 2 > 



z 2 5B 2 



slab 



2 Bl 



(19) 



This result if formally t 
small z values (see e.g. 



re same as the initial free-streaming result which can be found for 



1* 
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IV. FLRW FOR PURE 2D TURBULENCE 



In this Section, we shall follow the nonlinear formalism for FLRW proposed by Shalchi 
& Kourakis 



181] ). According to the results therein we have for pure 2D turbulence 

dk ± g 2D (k ± )e-H( A ^) k l. (20) 



dz 2 



2tt 
B 



o jo 



With the spectrum of Eq. (jSJ) we find 

;(Ax) 2 > 



dz 2 



j ? i-g $B 2D 

u i l 2D r>2 
n 



+ 



(21) 



The detailed calculation (limited to the case q = 0) was carried out in the latter reference, 
where a super diffusive behavior of the form <((Ax) 2 ) ~ z 4//3 was obtained. Our purpose in 
the following is to extend that result, for a general form of the wave spectrum. 

It can easily be demonstrated that the inertial-range of the spectrum yields a negligi- 
ble contribution in the rhs of (|21|) and was thus here neglected. The integral from the 
energy-range, extending the upper limit to infinity (l^p — > oo), can be expressed by Gamma 
functions 



dz 2 



;(Ax) 2 > 



diA-JBin /((Ax) 2 ) 



-I 



2D 



X 



B 2 
l-q 



(«-l)/2 



+ r 



q 1 



«(Ax) 2 )fcLn) : 



(22) 



Assuming that ((Ai) 2 )^ in 1 (i.e., the field-line MSD cannot exceed the maximum 
turbulence square length scale k^ n ), and using Eq. (fl3"j) we find 



dz 2 



;(Ax) 2 > 



4 2D r>2 
n 



l-q 



((Ax) 2 ) 



q-l 



(9-l)/2 



((Ax) 2 ) k 



2 

rain 



(l-9)/2' 



(23) 



The formula can be applied so long as < q < 2, except for q = 1. In the latter case, 
Eq. (I2T!) can be directly evaluated and we find a logarithmic behavior of the MSD. In the 
following, we shall further simplify Eq. fl23l . separately considering the cases < q < 1 and 
for 1 < q < 2. The relation e' = ((Ax) 2 ) k^^ <C 1 is assumed to hold everywhere. 
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A. The case < q < 1 



In this case the first term in Eq. (|23|) is dominant 



By making the ansatz ((Ax) y = a we can solve this ODE analytically. It can easily be 
demonstrated that 

b = ^~ (25) 
3 — q 



and 

2/(3-9) 



(26) 



Obviously we find 

\ < b < 2 (27) 
which is interpreted as superdiffusion. A diffusive behavior (b = 1) cannot be obtained. 



Interestingly, even for q = 0, one finds 6 = 4/3 (see in 181 ]) 



B. The case 1 < q < 2 

In this case the second term within brackets in Eq. (1231) is dominant (of the order 
e ( l -i)l 2 ^> i; see above) and we have 

^ {{Axf) ^J±J^ (kDkm ^. (28) 

By using Table [TT] for d% this can be simplified to 
and we finally find 

<(Ax) 2 > = (30) 
which is again the initial free-streaming (parabolic MSD) result. 



V. FLRW FOR SLAB/2D COMPOSITE GEOMETRY 

According to cosmic observations, it is more realistic than plainly adopting a pure-slab 



or pure-2D model, to consider a 20% slab/80% 2D composite model (see e.g. [19|). In this 



case, one rigorously obtains a 2nd-order ODE [cf. ([TO]) . (|20l) ]. whose RHS is the sum of the 
slab and 2D contributions, the relative weight of which is determined by the corresponding 
turbulence strength, i.e., 5B 2 slab /5B 2 and 5B 2 D /6B 2 . 

We shall now attempt to evaluate the asymptotic behavior of the field-line MSD in this 
hybrid (composite) model. 



A. The case < q < 1 

In this case we can combine Eqs. ( 1151) and (EMI) into: 

^<(A,) 2 > = 4c^-^r(l- g) sin(^ 

A _ q 5Bl D (((Ax) 2 )\^ l), \(l-q 

where negligible contributions were omitted in the rhs. Obviously this equation has the form 



+ lWS ™ r V1 , (31) 



2 



^((Axf)=az^ + P[((Axf)} { «- 1)/2 . (32) 

By applying the ansatz ((Ax) 2 ) = az b we find 

ab(b - l)z b - 2 = az^ 1 + (3a^/ 2 z b ^ 2 (33) 

where the definitions of a and (3 are obvious. It is straightforward to prove that, since b < 2, 
the second term in the rhs is dominant for z -— > oo. The slab contribution can therefore be 
neglected, so we can use Eqs. (!25j) and (|26l) also within the two-component model. 



B. The case 1 < q < 2 



In this case we can simply add the two contributions (Eqs. f|T9l) and fl3"0"]T 



v2 N _ 

1 2 5 n 2 



= -^2- (34) 



where we have set 

5B 2 = 5B 2 lab + 5B 2 2D . (35) 
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Geometry 


spectral index 


a 


b 


slab 
slab 
2D 
2D 

slab/2D 
slab/2D 


< q < 1 
1< q < 2 
0<q<l 
Kq<2 
0<q<l 
Kq<2 


q( q +l) l slab B g 1 ^ ^' V 2 J 
1 

k;i-^2M/ 2r (V)l 2/M 

D o V / J 

1 

2 B 2 

[d 1 z^^2a-^r(V)] 2/(3_9) 

1 <5B 2 

2 BJT 


g + 1 
2 

4/(3 - q) 
2 

4/(3 - q) 
2 



TABLE III: In this table, the results obtained for the parameters a and b, having adopted the form 
^(Ax) 2 ^ = az b for the field-line mean square deviation, are presented. In all (but one) cases, we 
find either super diffusion (1 < q < 2) or free-streaming (q = 2) of the field-lines. Diffusion (q = 1) 
can only be found for slab geometry and q = 0. 



VI. SUMMARY AND CONCLUSION 



We have investigated the random walk of magnetic field-lines for a more general spectrum, 
than the one employed in previous works. By exploring pure slab, pure 2D, and two- 
component turbulence models, we have calculated the field-line mean square deviation by 
applying the analytical description for FLRW proposed by Shalchi & Kourakis (see [ISl]). A 
super diffusive behaviour is found in all cases considered. In Table IHII the results obtained 
in this article are summarized. The only case where one obtains diffusion is for pure slab 
geometry and q = 0. As shown in this article the energy-range spectral index is a key-input 
parameter if FLRW is described. 

In the two-component turbulence model, which has been considered as a realistic model 
for solar wind turbulence (see [3]), we already find a weakly superdiffusive behavior if 
q = 0. For larger values of q we have ((Ax) 2 ) ~ ^ 4 /( 3 -'?). If the energy-range spectral index 
exceeds unity we find the same solution as in the initial free-streaming regime. Obviously, 
the energy-range spectral index has a very strong influence on FLRW behavior. 

From a theoretical point of view the results for pure slab geometry deduced in Section 3 
are very interesting and important because of two reasons: 

• for pure slab turbulence the parameter < (Ax) 2 > can be calculated exactly. No 
theory nor any ad hoc assumption have to be applied. 
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In all cases except q = we find superdiffusion of FLRW. 



Since in reality 20 % of the fluctuations can be represented by slab modes (see [19j) it is 
self-evident to assume that superdiffusion and not (classical or Markovian) diffusion is the 
regular case in astrophysical turbulence. 

If we merge from pure slab geometry to the slab/2D composite model a (nonlinear) theory 
has to be applied and an exact description of FLRW in no longer possible. By applying the 
ODE deduced by Shalchi & Kourakis (see 18J) we have shown that the superdiffusivity 
becomes even stronger in comparison to the pure slab results. 

It must be the subject of future work to apply these new results on realistic systems, 
such as solar wind turbulence. An important example is perpendicular transport of charged 
cosmic rays which is directly controlled by the FLRW, since charged particles are tied to 
magnetic field-lines (see flo)]). 
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